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Dedicated to the memory of I. M. Chakravarti 

Resolvable designs with two blocks per replicate are studied from 
an optimality perspective. Because in practice the number of repli- 
cates is typically less than the number of treatments, arguments can 
be based on the dual of the information matrix and consequently 
given in terms of block concurrences. Equalizing block concurrences 
for given block sizes is often, but not always, the best strategy. Suffi- 
cient conditions are established for various strong optimalities and a 
detailed study of E-optimality is offered, including a characterization 
of the E-optimal class. Optimal designs are found to correspond to 
balanced arrays and an affine-like generalization. 

1. Introduction. Block designs arise in comparative experimentation as 
fundamental devices for improving efficiency when working with heteroge- 
neous experimental units. The blocks are simply a partition of the units into 
(say) b sets exhibiting homogeneity within sets. Given the blocks, of sizes 
fci, fo, . . . , kf,, a block design is an assignment of v treatments to the Y^!j=i kj 
units. Optimality theory for block designs attempts to determine which of 
the many possible assignments is in some sense best. 

In some applications there are restrictions on the collection of possible as- 
signments. A block design is resolvable if the blocks can be partitioned into 
replicates, defined as sets of blocks with the property that each treatment 
is assigned to one unit in each set. The practical impact of, and motivation 
for, resolvability is to gain orthogonality between treatments and nuisance 
factors of concern. For instance, resolvability in sequential experimentation, 
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with replicates corresponding to time periods, is used to mitigate time ef- 
fects. Resolvability can likewise be useful in multi-site experiments and in 
experiments with multiple individuals handling experimental runs. Notably, 
the United Kingdom has for some time required the use of resolvable designs 
in agricultural field trials (see [20]). 

The combinatorial study of resolvability in block designs goes back at 
least as far as the well-known Kirkman's [14] schoolgirl problem. The notion 
entered the statistical lexicon with Yates' work on square lattice designs [29, 
30], although the term "resolvable design" was introduced by Bose [4]. Yates' 
lattice designs were extended to rectangular lattices by Harshbarger [11, 12]; 
see also [2]. Williams [27] and Patterson and Williams [21] introduced a 
large family of resolvable designs which they termed a-designs. Williams, 
Patterson and John [28] derived resolvable designs with two replicates from 
BIBDs (balanced incomplete block designs). Bailey, Monod and Morgan 
[1] proved strong optimality for the affine resolvable designs introduced by 
Bose [4]. Resolvable BIBDs have received significant attention in both the 
combinatorial and the statistical literature; for a summary, see Morgan [17], 
who surveys the major classes of resolvable designs with many references. 

With a-designs, Patterson and Williams [21] provided a flexible method 
for obtaining reasonably efficient resolvable designs for a wide range of values 
v, constant block size k and replication r. They also adapted their method 
to obtain resolvable designs with two different block sizes, k and k — 1, a 
first attempt at addressing the obvious restriction that resolvability with 
equal block sizes can be achieved only for v a multiple of k. John, Russell, 
Williams and Whitaker [13], in revisiting that idea, concluded that the a- 
technique for two block sizes could produce relatively inefficient designs for 
small v and recommended an interchange algorithm for construction of de- 
signs with better efficiency. John et al. [13] also discussed the practical need 
for resolvable designs with unequal block sizes; for example, about half of 
245 experiments examined by Patterson and Hunter [19] had unequal block 
sizes. See also [20]. 

To the authors' knowledge, the literature contains no systematic work on 
determining optimal resolvable block designs when block sizes need not be 
equal. This paper will undertake such work, for the special case of two blocks 
per replicate. Only with two blocks in each replicate must a block be large in 
the sense of containing at least half of the treatments. And with two blocks 
per replicate, the block sizes must be unequal for any odd v (although this 
work is not restricted to odd v). 

Let V(v,r;ki,k2) denote the class of all resolvable block designs with r 
replicates of v treatments, each replicate consisting of two blocks, of sizes 
k\ and k2, where k\ + ki = v. Of special interest is k\ = k, k<i = k — 1 for 
odd v = 2k — 1, but here no restrictions are placed on the two block sizes for 
the general theory. Let k\ denote the larger block size, k\ >k2, so that with 
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Table 1 
A resolvable design in 2?(9, 5; 5, 4) 



16 12 14 13 12 

27 43 25 24 35 

38 56 36 56 46 

49 78 78 87 79 

5 9 9 9 8 



no loss of generality | < k\ < v — 2. In the most common applications of 
resolvable designs, the number of treatments is large relative to the number 
of replicates; here, r < v — 1 is required, allowing optimality problems to 
be more easily attacked using the dual of the information matrix. These 
considerations frame the goal of this paper: to determine the best design 
d S T>(v, r; /ci, A^). An example of a resolvable design in T>(9, 5; 5, 4) is shown 
in Table 1, with the blocks written as columns. Later, this design will be 
proven optimal with respect to many useful criteria. 

For any resolvable design d£T>, ignoring the replicate grouping leaves an 
underlying simple block design for v treatments in 2r blocks. If the roles of 
blocks and treatments are reversed in this underlying design, another simple 
block design with 2r treatments in v blocks of size r is produced, called the 
dual design. 

Section 2 provides optimality background and establishes equivalences 
with the dual problem. Section 3 finds conditions for global optimality of 
designs having equality of block concurrences; for those who wish to jump 
ahead, the main results there are Theorems 4-7. Section 4 characterizes 
the E-optimal designs (Theorem 14) and finds the Schur-best of the E- 
optimal designs (Theorem 15). The main results are applied to special cases 
with h\ — k2 < 2 in Section 5 and designs for these cases are constructed in 
Section 6. Concluding comments appear in Section 7. 

2. Model, information and optimality criteria. Let y^ji denote the yield 
from the Ith. experimental unit in block j of replicate h. Thus, the triples 
(h, j, I) identify the experimental units and the design d corresponds to a map 
d[h,j, I] from the units to the set of treatments. The standard linear model for 
the yields incorporates a mean effect //, replicate effects ph, block effects fthj, 
treatment effects n and mean zero, uncorrelated, equivariable random 
error terms ehjl'- 

Uhjl = M + Ph + Phj +T~d[h,j,l] + e hjl, 

h = 1, . . . , r; j = 1,2; I = 1, . . . , kj. This model may be written in matrix 
terms as 



y = (il + (I r 8) l v )p + L/3 + X d r + e, 
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where, with the yhji lexicographically ordered in the vr x 1 yield vector y, 
the block incidence matrix is L vrX 2r = Ir ® Cn 1 ? fel ) an d the design matrix 

is the vr x v incidence matrix Xd, for which row (h,j,l) has a 1 in column 
i if and only if d[h,j,l] = i [i.e., unit (h,j,l) receives treatment i] and all 
other entries are zero. Replicate effects, block effects, treatment effects and 
error vectors are p rx i, /?2rxi> t vx \ and e vrX i, respectively. Choice of design 
is equivalently choice of X d - 

Linear model theory says that the information matrix for estimation 
of the treatment effects r is 

(1) C d = X' d [I - L{ULY l L'\X d = rl- NdDj 1 ^, 

where D s = L'L = Diag(A;i, k 2 , k\, k 2 , •■• , k\, k 2 ) is the diagonal matrix of 
block sizes. The v x 2r matrix Nd is the treament/block incidence matrix. 
The general off-diagonal element (NdN d )i^i of N d N' d is the number of blocks 
to which both treatments i and i' are assigned, called a treatment concur- 
rence count. Notice that the replicate incidence matrix {I r (8) 1„) plays no 
role in (1); the same form of information matrix is obtained for any sim- 
ple block design. It follows immediately that the information matrix for the 
corresponding dual design is 

(2) C dual = D S - -N' d N d . 

r 

The off-diagonal elements of N' d N d are block concurrence counts. 

All treatment contrasts are estimable with design d if and only if Cd has 
rank v — 1. Any such d is said to be connected; only connected designs are 
considered here. Most (but not all) commonly employed optimality criteria, 
including those to be used here, are functions of the v — 1 nonzero eigenvalues 
of C d . These will be ordered and labeled z d \ < z d2 < • • • < Zd jV -i- 

Much of the optimality work below focuses on minimizing functions of 
the form 

v-l 

(3) ipf(zd) = ^2f(zdi), 

i=X 

where / is convex and Zd is the vector of nonzero eigenvalues. If / in (3) is 
Schur-convex (Schur-convex functions include the convex functions; see [3], 
Section II. 3), then (3) is said to be a Schur- criterion. Design d\ is Schur- 
better than design d 2 (d 2 is Schur-inferior to d\) if ipf(z dl ) < ipf(zd 2 ) for 
all Schur-convex / with strict inequality for at least one such /. A design 
optimal with respect to all (i.e., minimizing all) Schur-convex criteria is said 
to be Schur- optimal. If / in (3) satisfies (i) / is continuously differentiable 
on (0, maxdev tr (Cd)) with /' < 0, /" > and < and (ii) lim^o f(%) = 
00, then (3) is said to be a type-1 criterion (see [7]). A design optimal 
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with respect to all type-1 criteria is said to be type-1- optimal. One popular 
criterion belonging to both families just defined is the ^4-criterion specified 
by fi x ) = ■ A criterion not of the form (3) (though it can be written as a 
limit of such criteria) is 



(4) 



ll>E(z d ) = , 

Zdl 



called the E-criterion. An E-optimal design minimizes (4) or, equivalently, 
maximizes Zdi ■ For a broader discussion of optimality criteria and their sta- 
tistical meanings, see [24]. 

For the current endeavor, it is advantageous to approach the z<n, and con- 
sequently (3) and (4), through Cd ua i- How this is done is shown in Lemma 1 
below. Let <pdhh' be the block concurrence for (i.e., the number of treat- 
ments common to) the blocks of size k\ in replicates h and h! . For h ^ h' , 

k 2 

^ ^*dhh' = 'Pdhh' — ki+k2 anc ^ define the (symmetric) optimality matrix 

by 



/ kik 2 



(5) 



M d 



Vdl2 

k\k 2 
v 



V 



#03 • 


■ <P*dlr 








■ #Br 


k\k 2 




V 


■ <P*dZr 




k\k 2 



with eigenvalues e^i > eau > • • • > edr ■ 



Lemma 1. For any d GBfi)^;^,^), the eigenvalues of Cd areO, max{0, 
v - r - 1} copies of r and (r - %^,r - ^g, . . . ,r - jj^). Consequently, 
a ipf -optimal design minimizing (3) will equivalently minimize Y7h=ifi r ~ 
|^) and an E-optimal design minimizing (4) will equivalently minimize 
edi- 



Lemma 1 is proved in Appendix A.l. Provided r < v — 1 (as earlier re- 
quired in Section 1), Cd has v — r — 1 eigenvalues fixed at r. Working with the 
dual through Md not only makes this evident, but allows these structurally 
fixed eigenvalues to be easily set aside. 

With the optimality problem recast in terms of and its eigenvalues, a 
crucial concept for proofs of Schur-optimality is now defined. 



G 



J. P. MORGAN AND B. H. RECK 



Definition. Let {xj}f =1 and {yi}f =1 be nonincreasing sequences of real 
numbers such that Ya=i x i = J27=i Hi- ^ 
i i 
^ Xi < yi for all 1 < I < n 
i=i i=i 



or, equivalently, 



n-l+l n-l+l 

^ Xi > V] yj for all 1 < I < n, 



then {yi}f = i is said to majorize {xj}™ =1 . 



The importance of majorization is evident in the following result (see, 
e.g., [3], page 40). 

Theorem 2. Let {xj}™ =1 and {?/j}™ =1 be nonincreasing sequences of real 
numbers such that Y2=i x % = E?=i Hi- Then Ef=i f( x i) < E"=i fiVi) for all 
real-valued convex functions f if and only if {?/j}™ =1 majorizes {xi}f =1 . 

If the sequences {xj}f =1 and {?/i}" =1 are written as the elements of vectors 
x and y, then the fact that y majorizes x is written as y >~ x or x -< y. 

Let e,i and e^* be the vectors of eigenvalues for the optimality matrices 
for designs <i and d* , respectively. Lemma 3 lists simple majorization facts 
([3], page 30) used in subsequent sections. For majorization comparisons, 
the third of these allows work directly with the edh rather than the r — 
(see Lemma 1). Thus, if ed >- e<i* and the two vectors are not identical, then 
d* is Schur-better than d. Design d* is Schur-optimal if a >- ed* for every 
deV. 



Lemma 3. For real numbers {xj}™ =1 and {yi}™ =1 with Y^i=i x i = Ya=i Vi> 



(i) if xi > x 2 = x 3 = • • • = x n and y\>x\, then {yi}™ =1 majorizes {xj}f =1 ; 

(ii) if x\ = X2 = • • • = x n _\ > x n and x n >y n , then {yj}™ =1 majorizes 

{ x i}f=l/ 

(hi) i/{yi}?=i majorizes {xj}f =1 , then {a — ^}?=i majorizes {a — ^f}f =1 
for any real a, b. 



3. Equal concurrence designs and global optimality. Among simple block 
designs, the BIBDs are Schur-optimal, a result which follows from equality of 
treatment concurrences inducing complete symmetry of the information ma- 
trix. The analogous notion for duals is equality of block concurrences, which 
this section explores for utility with resolvable designs. A resolvable design 
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d 6 T>(v, r; k%, ^2) having block concurrence counts 4>di2 — </>(fl3 — ^d23 = • • • = 
4>d,r-i,r = $ f° r some fci — ki < 9 < k\ is called an equal concurrence design 
with common concurrence 9, or ECD(#). The subsections below will explore, 
in turn, (i) ECDs with equality of eigenvalues in the optimality matrix, (ii) 
other ECDs which can be proven Schur-optimal and (iii) Schur-domination 
of designs with unequal block concurrences by one or more ECDs. 

3.1. Schur- optimality via equality of eigenvalues. For an ECD(#), the 
optimality matrix (5) is 

\P fa k l 



(6) M d 



v 



where / is the r x r identity matrix, J is the r x r matrix of ones and 
p = k\k2 is the product of the block sizes. Like a BIBD information matrix, 
it is completely symmetric, but unlike that matrix, M d for an ECD(#) is 
nonsingular and thus can have two distinct, relevant eigenvalues, rather than 
just one. 

Theorem 4. Suppose that T>(v,r;ki,k2) is a resolvable design setting 
for which (k± + k2)\kf and define 

k 2 k 2 
k\ + k2 v 

Then ECD(0* )' ' s in V are Schur-optimal whenever they exist. 

k 2 

Proof. Given the conditions on T>, the inequalities k\ — k2< ki+k 2 — ^ 1 
imply that 9 = 9* is an admissible value for the common block concurrence 
of an ECD(6>). For 9 = 9*, M d in (6) is 2 J. Since the eigenvalues of M d are all 
identical, they are majorized by the eigenvalues of every competing design. 
□ 



Corollary 3.4 of [1] established that affine-resolvable designs are Schur- 
optimal. Theorem 4 generalizes that result when there are two blocks per 
replicate. Here, the optimality condition is that all pairs of blocks of size 
k\ have the same concurrence (7). When k\ = k2 and 2\k\, ECD(#*)'s are 
affine-resolvable designs. It is obvious that (6) has two distinct eigenvalues 
for any 9^9*. 

Example. Consider the setting V (9, 4; 6,3). Then (ki + k 2 )\kl with 9* = 
4 and if an ECD(4) exists, it is Schur-optimal. In fact, an ECD(4) does exist 
and is shown in Table 2. 

The settings for which k\ is a multiple of k\ + /C2 are relatively sparse 
(a situation much like that of BIBDs relative to all simple block design 
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Table 2 
A Schur-optimal ECD(4) 
in £>(9,4;6,3) 



1 


7 


1 


5 


1 


3 


3 


2 


8 


2 


6 


2 


4 


4 


3 


9 


3 


9 


5 


9 


5 


4 




4 




6 




6 


5 




7 




7 




7 


6 




8 




8 




8 



settings). For the 1225 pairs 2 < k 2 < k\ < 51, only 23 meet the divisibil- 
ity requirement implied by (7). Theorem 4 is thus only a start, albeit an 
important one. 

3.2. Global optimality of other ECDs. Good designs are expected to have 
eigenvalue structures "close" to that of ECD(6 l *)'s, suggesting this question: 
Is some equal concurrence design Schur-optimal when {k\ + k2)fkf? To in- 
vestigate this question, define the block concurrence parameter 9 by 

/ k 2 

(8) £ = in t(_^ 
and write 

(9) 7 



h + k 2 
k 2 



h + k 2 

Then < 7 < 1 and a necessary condition for existence of ECD(#*) is 7 = 0. 
Consequently, 7 is called the block discordancy coefficient; it measures the 
departure of the block sizes from that required for equality of all eigenvalues. 
The parameter 7 will play a pivotal role in the remainder of this paper, as 
will the concurrence discrepancies defined next. 

Define the block concurrence discrepancy matrix = (5dhh'), where 



Jdhh' 



4>dhh> -0, if h ^ h', 
0, if h = ti. 



For h! , the off-diagonal elements b^hh' will be referred to as block con- 
currence discrepancies. Rewritten in terms of block discrepancies and the 
discordancy coefficient, the general optimality matrix (5) becomes 

(10) M d = ^J- 7 (J-/) + A (i . 

v 

For ECD(0*), 7 = 0, A d = and M d = gi". If 7 / 0, (10) shows that the form 
of optimal design may depend both on the magnitude of 7 and the values 
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of the discrepancies. Matrix (10) is completely symmetric for, and only for, 
ECDs, in which case A^ is an integer multiple of J — I. 

The ECDs which are combinatorially closest to ECD(#*) are those with 
either 9 = 9 or 9 = 9 + 1. Not surprisingly, these are strong competitors in the 
optimality race. ECD(#)s have 4>dhh' = 9 for each l<h^h'<r, so A^ = 
and the eigenvalues of are 

(11) 6 ( 7 ) = ^ + 7 an d &( 7 ) = P - ( r - 1) 7 , 

v v 

with frequencies r — 1 and 1, respectively. The two eigenvalues satisfy £i (7) > 
6(7)- If <t>dhh' = 9 + 1 for every h / ti, then ECD(# + l)'s have A d = (J - I) 
and the eigenvalues of are 

6(7- 1) = - - (1 -7) and 

(12) 

6(7"l) = ^ + (r-l)(l-7), 
with frequencies r — 1 and 1, respectively, and with £2(7 — 1) > £1(7 — !)• 

Theorem 5. For < 7 < \, ECD(9 )'s are type-l- optimal and E-optimal. 

Proof. The eigenvalues of the information matrix for any design d € 
V(v,r;ki,k2) are < z&\ < Zdz < • ■ • < z^r and v — r — 1 copies of r, and 
Si=i -Srfi = r ( r ~~ 1) i s constant for all designs in T>. For an ECD(#), call it d, 
the Zfc, following from (11), have the form z dl = z^ 2 = ■•• = zj v _ 2 < zg v _ x . 
Theorem 2.3 of [7] thus gives the result if it can also be shown that d 
minimizes £)i=i z 2 di over V. 

For d G T> (v, r; k\, fo) with optimality matrix {Mdjhh 1 = i&dhh' ~ l) having 
trace trM rf = 2f , 



2 _ / 1 ^ 2 , / ue <^ 



trQ 2 = E4 = (----l> 2 + E 



(t? - l)r 2 - — tr M d + — tr Mj 



„2 _ 



= (^-3)r 2 + r+^EE(^-7) 2 , 
p h<h' 

so that trCj is minimized by designs that minimize YjJ2h<h'(^dhh' — l) 2 ■ 
Since Sdhh' is integral, the unique minimum of tr Cj on < 7 < \ is at 
6 dhh t = 0, achieved only by ECD(#). For 7=31 anv values bdhh' £ {0,1} 
minimize trCj. □ 
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Now define the F-criterion as the value of the largest eigenvalue of Cd 
that is not constrained by the setting to equal r, that is, 

/ (n \ vedr 
WF{Cd) =z dr =r . 

P 

Minimizing ipF(Cd) ° ver 2? is equivalent to maximizing e^r- This criterion 
can be important in establishing Schur-optimality, as shown next. 

Theorem 6. An ECD(O) is Schur-better than a competitor with a dif- 
ferent set of eigenvalues if and only if it is F-equivalent or better than that 
competitor. Consequently, ECD(9)'s are Schur-optimal if and only if they 
are F- optimal. 

Proof. Follows from (11) and Lemma 3(i). □ 

A result of a similar flavor holds for ECD(# + l)'s using the E-criterion. 
As pointed out by Kunert [15], page 385, designs with eigenvalues Zdi in the 
form of Lemma 3(ii) are Schur-best whenever they are E-optimal. For the 
current problem, this is stated as follows. 

Theorem 7. An ECD(9 + 1) is Schur-better than a competitor with a 
different set of eigenvalues if and only if it is E-equivalent or better than 
that competitor. Consequently, ECD(9 + l)'s are Schur-optimal if and only 
if they are E-optimal. 

In establishing necessary and sufficient conditions for Schur-optimality 
for their respective ECDs in terms of a single eigenvalue, Theorems 6 and 7 
provide simple tests for comparing these designs to any other design. An im- 
mediate consequence is that an ECD(#) with 9 ^ {9, 9 + 1} is Schur-inferior 
to at least one of these two competitors. Thus, among ECDs, only these 
two competitors remain. They are compared to one another in the following 
corollary. 

Corollary 8. ECD(9)'s are Schur-better than ECD(9 + l)'s if and 
only if j < i, and ECD(9 + l)'s are Schur-better than ECD(9)'s if and only 

ifl> r -^- 

Proof. Simply use (11) and (12) in applying Theorems 6 and 7. □ 

Though design nonexistence can play a role, Corollary 8 says that one 
does not expect to find a Schur-optimal design for £ < 7 < (Theorem 5 
makes for an interesting juxtaposition). Schur-domination can nonetheless 
be used to eliminate many competitors, as shown next. 
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3.3. Schur-inferiority of designs lacking equality of block concurrences. 
Given the results of the preceding subsection, the remaining question from a 
global optimality perspective is if (and when) designs outside the ECD class 
can be preferable. This question can be effectively pursued by application 
of Theorems 6 and 7, once workable bounds for the largest and smallest 
eigenvalues e d \ and e dr of the optimality matrix are in place. 

Lemma 9. Let d G T>(v,r;ki,k 2 ) have concurrence discrepancy matrix 
Ad = (o~dhh') a nd optimality matrix M d . Then: 

(i) $di2 < e dl > f + 7 - 5 dl2 and e dr < § - 7 + 5 dl2 ; 

(ii) 6 dl2 > e dl > £ - 7 + 5 dl2 and e dr < 2 + 7 - 5 d i 2 . 

Proof. The leading 2x2 minor of M d , which is M dll = (2 + 7 — 5 d i 2 )I- 
(7 - 5 d i 2 )J, has eigenvalues f + 7 — 8 d i2 and ^ - 7 + <5di2- A Sturmian sep- 
aration theorem ([22], page 64) provides the bounds. □ 

Define design d to be an almost equal concurrence design, or AECD, if 
4>dhh' £ {0)0 + 1} for a h h^h' and if each value is attained for some h,h! . 
If any <p d hh' is not in {9, 9 + 1}, then d is an unequal concurrence design, or 
UCD. In terms of discrepancies, AECDs have all 6 d hh' £ {0> 1}> while UCDs 
have some 5dhh> < — 1 ° r > 2. Depending on 7, AECDs can be optimal in at 
least some senses (see Section 4), necessarily ruling out global optimality of 
ECDs for some 7. The next few results will show that ECDs often dominate 
UCDs. 

Corollary 10. Suppose d G T>{v, b; k\,k 2 ) is a UCD with 5 d hh' < -a 
for some 1 < h 7^ h! < r and some integer a > 1. Then: 

(i) ECD(9)'s are Schur-better than d if 7 < ; 

(ii) ECD(9 + 1 /a are Sc/rar better to d i/ 7 > 

Proof. For the UCD d as described in the corollary, take 5 d i 2 < —a. 
Then from Lemma 9, e d \ > 2 + 7 — « and jj — 7 + a > e dr . By Theorem 6, 
an ECD(0) is Schur-better than d if 6(7) > £ - 7 + « > e dr <=^ 7 < 
By Theorem 7, an ECD(# + 1) is Schur-better than d if e dX > 2 + 7 - a > 
6(7-1) ^7> r ^ 1 - □ 

COROLLARY 11. When r < 4, aZZ UCDs with some 5 d hh' < — 1 are Schur- 
inferior to an ECD, and when r = 5 or 6, UCDs with some 5 d hh' < — 2 are 
Schur- inferior to an ECD. 

The next two corollaries are similarly shown. 
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Corollary 12. Suppose d e V(v,b;ki,k2) is a UCD with 5dhh' > « f or 
some 1 < h ^ h' < r and some integer a>2. Then: 

(i) ECD(0ys are Schur-better than d if 7 < j; 

(ii) ECD (6 + l) , s are Schur better than d if 7 > r ^ i . 

Corollary 13. When r < 4, a// {/CDs with some bdhh' > 2 are Schur- 
inferior to an ECD, and when r = 5 or 6, UCDs with some b~dhh' > 3 are 
Schur-inferior to an ECD. 

Corollaries 11 and 13 say that optimal (with respect to any convex cri- 
terion) designs in settings T>(v,r;ki,k2) with r < 4 must be an ECD(#), an 
ECD(# + 1) or an AECD. Optimal designs in settings with r = 5 or 6 must 
have block concurrence discrepancies ddhh' G { — 1,0, 1, 2} for all 1 < h 7^ h! < 
r. It is unlikely that such global statements can be much improved. The im- 
portance of these results lies in the prevalence of small r in the application 
of resolvable designs. 



4. E-optimality. Sufficient conditions, which are necessary given exis- 
tence, will be developed for E-optimality of designs in T>(v,r;ki,k2)- The 
main results, Theorems 14 and 15, will be stated after first defining a sub- 
class where E-optimal designs will be found. 

Corollary 16 below will remove the UCDs from E-contention, so that only 
ECD(^)'s, ECD(0 + l)'s and AECDs need be considered. These designs have 
all 

fidhh' ^ {0, 1} and so have which is the adjacency matrix of a simple, 
undirected graph on r vertices. Any of these designs for which — J is 
(with suitable ordering of replicates) of the form 

f-Jh \ 

(13) 

V -J J 

for positive integers n and t\<t%< - ■ • <t n U = r ) i s sa id to be group- 

affine. For group-affine designs, concurrences cpdhh' are constant (= 9) within 
groups of sizes t±, . . . ,t n and are constant (=6 + 1) between groups. A group- 
affine design is said to be uniform if, for given n, the range of group sizes 
ti is at most one. For any group-affine design with the number of groups 
n < riry = int(j^-), let t^ d ) denote the vector of its group sizes ti arranged 
in increasing order and with ra 7 elements, padding with zeros as necessary. 
For example, if r = 7, 7 = | and d has ti values 1, 3 and 3, then n 7 = 4 and 
t{d) = (0,1,3,3). Now the main results can be stated. The phrase "E-Schur- 
optimal" means "Schur-optimal within the class of all E-optimal designs." 
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Theorem 14. Any group-affine design is E-optimal if n<n^ and 7 < 
J —p^- If any such design exists, then all such designs form the class V>e of 
all E-optimal designs. ECD(9 + 1 )'s are Schur-optimal for 7 > . 

Theorem 15. For d\,d<2 E T>e, if tu^ -< t/ d2 )i then d\ is Schur-better 
than di- Design d* 6 T>e is E- Schur-optimal if it is uniform with n = n 7 . 

Theorem 14 characterizes the class of E-optimal block concurrence struc- 
tures for T>(v,r;ki,k2). Although not much discussed in the literature, the 
class of E-optimal designs in a given simple block design setting often con- 
tains a variety of designs with different information matrices (see, e.g., [18] 
and [26]). Theorem 14 reveals the same situation for resolvable design set- 
tings. This allows other design criteria to be brought to bear: one should 
choose the best of the E-optimal designs. Theorem 15 does this. In a very 
strong sense, E-Schur-optimal designs are the best of the E-optimal designs. 

The proofs of Theorems 14 and 15 depend on a series of technical results 
that are developed in the remainder of this section. So as not to overly 
disrupt the flow of the main line of proof, the longer of the "sub-proofs" are 
delayed until Appendix A. The first task is to rule out UCDs. 

Corollary 16. For all r > 2 and < 7 < 1, ECD(6)'s are E-better 
than UCDs. 

Proof. The maximum eigenvalue of the optimality matrix for an ECD(#) 
is £1(7) = fj + 7- For the UCD d, suppose that 6^12 < —ot for some integer 
a > 1. Then by Lemma 9, e dl > f + 7 - <S dl2 > £1(7) and ECD(#)'s are E- 
better than d. If ^12 > ct for some integer a > 2, then another application of 
the lemma gives e<ji > jj —7 + ^12 > £1(7) an d, again, ECD((9)'s are E-better 
than d. □ 

Completing the proof of Theorem 14 is a matter of minimizing e^i over 
all d for which every 5dhh' £ {0,1}- Before providing the details, here is a 
sketch of what will be done. First, a lower bound for e^i is found in terms 
of the largest eigenvalue of a specific projection PA^P of the discrepancy 
matrix. This bound implies that ECD(^)'s are E-superior to all designs 
for which PA^P has a positive eigenvalue. The next step is thus to find 
a necessary and sufficient condition on A^ so that PA^P is nonpositive 
definite; this condition turns out to be exactly the definition of group-affine 
design. Finally, the general form of for group-affine designs is examined 
in detail to determine which of these competitors are E-optimal, producing 
the conditions of Theorem 14. 

Lemma 17. Ford e T>{v,r;k\,k2) and with P = (I —^J), letudi andu^r 
be the maximum and minimum eigenvalues of PA^P, respectively. Then: 
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(i) if Udi > then e d i > % + 7 + u d\ ; 

(ii) e dr < 2 +J + u dr . 

The proof of Lemma 17 appears in Appendix A. 2. Lemma 17 combined 
with Theorems 6 and 7 immediately gives the following corollary: 

Corollary 18. For de V(v,r;ki,k2) and with P = (I — ij), let u d \ 
and u dr be the maximum and minimum eigenvalues of PA d P. 

(i) If j < -^f-, then ECDfej's are Schur-better than d. 

(ii) // u dl > and 7 > { r ~ u f~ l ), then ECD(6 + l) , s are Schur-better 



(iii) If u d \ > 0, then ECD(9) J s are E-better, but not necessarily Schur- 
better, than d. 

Parts (i) and (ii) of Corollary 18 provide alternative routes (cf. Corollar- 
ies 10 and 12) for establishing Schur-domination of ECDs. Part (iii) is the 
key part of the proof of Theorem 14. If d is not to be eliminated by ECD(0), 
then PA d P cannot have a positive eigenvalue, that is, must be nonpositive 
definite. Recall that the current E-competitors are all d for which is the 
adjacency matrix of a simple, undirected graph. 

Lemma 19. Let A be the adjacency matrix for a simple undirected graph 
of r vertices. PAP is nonpositive definite if and only if A — J may be writ- 
ten (possibly after vertex permutation) in the form (13) for some positive 
integers n and t± < t2 < ■ • ■ < t n with Ya=i U =r - 

Combining Lemma 19 with Corollaries 16 and 18(iii), it has now been 
shown that existence of an ECD(#) implies that an E-optimal design is a 
group-affine design. Group-affine designs include the ECD(#)'s (put n = 1, 
A d = 0) and the ECD(# + 1) (put n = r, A d = J — I) at the extremes for n. 
The proof of Lemma 19 appears in Appendix A. 3. 

It remains to determine which of the group-affine designs are actually E- 
best, requiring knowledge of e d \ for this class. For any group-affine design, 
write M d = + 7)/ + A d - jj, where, for given t\ < • • • < t n , 



than d. 



A d - 1 J = { 



( l J ti,ti 
(7-l)^ 2 ,t 



{l-l)J tl M 

l J t 2 ,t 2 



(7-l)^A\ 

(7- l)^2,*n 



(14) 



l)H d . 



l J t n ,t 



J 
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An E-optimal design will maximize the minimum eigenvalue of Hd- Clearly, 
Hd has J2?=i(U — I) = r — n eigenvalues of zero (corresponding to eigen- 
vectors which are orthogonal contrasts within groups of sizes t\, . .. ,t n ). So 
all of these designs have Hd with at least one eigenvalue of zero, except 
ECD(# + 1), for which the eigenvalues of Hd are 1 (frequency r — 1) and 
1 + r(7 - 1). ECD(fl + 1) is therefore E-optimal if 1 + r(-y - 1) > 0, that is, 
if 7 > ^r- (in which case it is Schur-optimal; see Theorem 7). Otherwise, all 
designs for which Hd is nonnegative definite are E-optimal. Needed now are 
the eigenvalues of Hd other than zero. Let Dt be the diagonal matrix with 
diagonal elements (t±, . . . ,t n ). Lemma 20 is proved in Appendix A. 4. 

Lemma 20. The eigenvalues of Hd specified by (14) are (with fre- 

1/2 1/2 

quency r — n) and the eigenvalues of D t ED t , where E = I n — (1 — 7) J n . 

1/2 1/2 

Consequently, Hd is nonnegative definite if and only if D t ED t is non- 
negative definite. This is so if and only if E is nonnegative definite, that is, 
if and only if 1 — 7 < -. Thus, for any 7 < ^— the E-optimal designs are 
all group-affine designs for which the number of groups n is no larger than 
j4— , that is, for which n < n 7 . This completes the proof of Theorem 14. 

To prove Theorem 15, recall the definition of t(^) given just prior to Theo- 
rem 14 and note that now only d G T>e is being considered. It is shown in Ap- 
pendix A. 4 that the eigenvalues of Hd, aside from r — n 7 zeros, are the eigen- 
values of E l l 2 D t E 1 / 2 . Thus, the problem is to show that tui) ~< £(d 2 ) implies 
the eigenvalues of E 1 ! 2 D t{di) E l l 2 are majorized by those of E 1 l 2 D t{d2) E l l 2 . 

Now, tu{\ -< tfda) i s equivalent to tr^) = St(d 2 )i where S is a doubly 
stochastic matrix; S = o,%Qi for permutation matrices Qi and positive 
numbers a, summing to 1 ([3], pages 33 and 37). Thus, 



£1/2 



SotQiA (d2) Qi 



E 1 ' 2 



^iQAE^D^E^Qi 



the last equality following because E 1 ! 2 = (I — ^ J) + y i — (1 — 7) J com- 
mutes with any permutation matrix. This shows that E l l 2 Dt {d ^E l l 2 is a 
symmetrization of E l / 2 D t ^E 1 / 2 and thus the eigenvalues of the first ma- 
trix are majorized by those of the second (this follows from [3], page 69). 
Since t^d*) -< t^ d ) for every d S T>e, d* is E-Schur-optimal and Theorem 15 is 
proved. 
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5. Special cases: (fci — /C2) < 2. In this section, the three important 
special cases of k\ and k 2 being equal or nearly so, k 2 £ {k±, k\ — 1, k\ — 2}, 
are investigated in light of the results in Sections 3 and 4. Put k 2 = k\ — m 
so that (k\ — k 2 ) < 2 says m = 0, 1, or 2, and for any m, 

fc 2 fc 2 _ fci m m 2 

1 ' ki + k 2 ~ 2k l -m~ 2 4 4(2^ - m) ' 

A: 2 - 

Recall that 9 is the integer part of (15). The values for 7 = kl + k2 — 9 in the 
corollaries below are easily found using (15). 

Corollary 21. For ki = k 2 , 

(i) if2\k\, then 7 = 0, (A;i + k 2 )\k\ and ECD(9* )' 's are Schur- optimal; 

(ii) if2\Jk l; thenj = l and E CD (9)' s are E-Schur and type- 1- optimal. 

When k\ = k 2 and 2\k\, Schur-optimality also follows from [1], Corollary 
3.4. For 2[//ci, the result follows from Theorems 5 and 15. 

Corollary 22. For k\ — k 2 = \, 

(i) if2\k\, then 7 = and ECD(9)'s are E-Schur and type-1- optimal. 

(ii) if 2\jk\, then 7 = ^ and uniform group-affine designs with four 
groups are E-Schur- optimal if r > 5. For r < 4, ECD(9 + l)'s are Schur- 
optimal. 

As an example, the design in Table 1 is E-Schur-optimal for 9 treatments 
in 5 replicates with block sizes 4 and 5. The design consisting of the first 
four replicates is Schur-optimal. 

Corollary 23. For ki - k 2 = 2, 

(i) if2\k\, then 7 = ^r, and uniform group-affine designs are E-Schur- 
optimal if the number of groups is 3 for v = 6 and 2 for v > 10. For r = 2, 
ECD(9 + 1) 's are Schur-optimal; 

(ii) z/2|//ci, then 7 = and ECD(9)'s are Schur-optimal. 

The Schur-optimality in part (ii) of Corollary 23 follows from Theorem 6 
and Lemma 9. 

6. ECDs, balanced arrays and design construction. Balanced arrays were 
introduced by Chakravarti [5, 6] as a useful device for fractional factorial 
designs and they have since been investigated by a plethora of authors, 
including, of late, Kuriki [16], Fuji-Hara and Miyamoto [9], Sinha, Dhar 
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and Kageyama [25] and Ghosh and Teschmacher [10]. Here, only strength- 
two arrays on two symbols will be needed. A balanced array of strength 2, 
BA(iV, m,2), on the symbols and 1, is an N x m array with the property 
that for any selection of two columns, the JV pairs formed by the rows are 
(0,0), (0,1), (1,0) and (1,1), with frequencies [1q, /Ui, fii and [i 2 , respec- 
tively. Two-symbol orthogonal arrays of strength two are the special case 
Mo = Mi = A*2- When the number of 0's in each column is specified (as is 
the case below), the /jj's are all determined by 9 = /Jo and the array will be 
denoted BA(JV, m, 2; 0). 

Bailey, Monod and Morgan [1] demonstrate the combinatorial equivalence 
between affine resolvable designs and orthogonal arrays. Their method can 
be used to express any resolvable design as a combinatorial array, as follows. 
Given a resolvable design for v treatments in r replicates, each consisting of 
s blocks, label the blocks within a replicate 0, . . . , s — 1. Construct a v x r 
array by identifying rows of the array with treatments of the resolvable 
design, and columns with replicates: symbol j £ {0, 1, . . . , s — 1} is placed 
in row i, column h if and only if treatment i is in block j of replicate h. 
Theorem 24 is evident. 

Theorem 24. Each ECD(6 ) in V{v, r ; k\ , k 2 ) is equivalent to a BA (v,r,2;9). 

A group-affine design with n groups of sizes t\,...,t n is a juxtaposition 
(BAi, BA2, . . . , BA n ) of balanced arrays BAj = ~BA(v, ti, 2; 9) so that any 
two columns from different BAj's form a BA(w, 2, 2; 9 + 1). Call such a juxta- 
position a grouped balanced array, denoted by GBA(v, (t±, . . . , t n ), 2; 9). Con- 
structions for the designs in Corollaries 21-23 of Section 5 are now listed. 
These are stated in terms of Hadamard matrices, that is, orthogonal ma- 
trices for which every element is 1 or —1. The Hadamard conjecture says 
that a Hadamard matrix exists for every order a multiple of four. Existence 
is known for all such orders up to 664 and for infinitely many other values 
(see [8]). A Hadamard matrix is said to be standardized if the first row and 
column are all ones; this can always be achieved. 

Theorem 25. If k\ = k 2 is even, then a Schur- optimal ECD(9*) cor- 
responds to an OA(v,r, 2; ^). Existence of a Hadamard matrix of order v 
implies existence of the OA for every r < v — 1. 

Theorem 26. If hi = k 2 is odd, then a type-1- optimal ECD(9) corre- 
sponds to a BA(v, r, 2; ). Existence of a Hadamard matrix of order v + 2 
implies existence of the BA for every r < |. 

Proof. The value of 9 is int( ) = Given the standardized 
Hadamard matrix, permute columns (except the first) so that the second 
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row has 1 in the first columns. Delete the first two rows and the 
first ^rr^ columns, then replace —1 by throughout. Clearly, the result is 
BA(t/,§,2;^). □ 

Theorem 27. If k\ = ki + 1 is even, then a type-1- optimal ECD(9) 
corresponds to a BA(v,r,2; ^-)- Existence of a Hadamard matrix of order 
v + 1 implies existence of the BA for every r < v . 

Proof. The value of 6 is int( (( " +1 J /2)2 ) = ^±±. Given the standardized 
Hadamard matrix, delete the first row and column, then replace — 1 by 
throughout . The result is B A (v , v, 2; ^) . □ 

Theorem 28. // k\ = k 2 + 1 is odd, then a type-l- optimal ECD(8 + 1 ) 
with r < 4 corresponds to a BA(v,r,2; ^^-), which always exists. For r > 5 
and v > 9, an E-Schur- optimal group-affine design corresponds to a 
GBA(v, (ti, t 2 , t%, £4), 2; ) with £4 — t\ < 1. Existence of a Hadamard ma- 
trix of order v + 3 with a 4 x r submatrix of the form 



(I" 




1*,, 


"1*4 


\ 






"1*3 


1*4 


K 


-if. 


1*3 


1*4 




V-i 4l 


1*2 


1*3 


1*4 





implies existence of the GBA. 

Proof. The value of 9 is mt( i[v+1 ] /2)2 ) = ^f±. For r < 4, take OA(> - 
1,7", 2) on {0,1} and add one row of 0's to get BA(u,r,2;^). Given the 
assumed Hadamard matrix, delete v + 3 — r columns not containing the 
submatrix, then delete the submatrix and add a row of l's. Replacing —1 
by throughout gives the GBA. □ 

The maximum r admitted for a given v by the Hadamard construction in 
Theorem 28 depends on the particular Hadamard matrix chosen. All noniso- 
morphic Hadamard matrices are known up through order 28 and, for these, 
a complete search has produced (v, r) = (9, 5), (13, 5), (17, 9), (21, 9), (25, 13); 
the first four of these appear in Appendix B. A search of a few known 
Hadamard matrices (compiled by N. J. A. Sloane at www.research.att.com/~njas/hadamard/) 
of orders up to 48 has further produced (v, r) = (29, 13), (33, 16), (37, 18), (41, 24), (45, 19). 
Clearly there is room for more work to be done here. 

Theorem 29. // k\ = k 2 + 2 is even, for v > 10 an E-Schur- optimal 
group-affine design corresponds to a GBA(v, (ti,t2),2; ) with t2—t\< 1. 
Existence of a Hadamard matrix of order v-\-2 implies existence of the GBA 
for every t\ + t2 = r < 5 + 1. 
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Proof. The value of 9 is int( (( " +2)/2) ) = ^±2. Given the standardized 
Hadamard matrix of order v + 2 = Ah, permute columns (except the first) 
so that the first three rows are 

flh lh U U \ 

l ft -lft lft -1ft . 
\lft -lft -lft lft / 

Now delete the first 2/i columns, delete the first three rows, then add one row 

of — l's and finally replace —1 by throughout. The result is GBA(t>, (^-r, ^ip); 2; 

with two juxtaposed BA(u, ^^j 2 -, 2; being the first /i and last /i 
columns. For smaller r, delete int(/i — |) columns from the first component 
BA and int(/i — from the second. □ 

Theorem 30. If kx = k% + 2 is odd, a Schur- optimal ECD(9) corre- 
sponds to a BA(v,r,2; : ^-). Existence of a Hadamard matrix of order v 
implies existence of the BA for every r < v — 1. 

PROOF. The value of 9 is int( (( " +2)/2)2 ) = 2±4. Given the standardized 
Hadamard matrix, delete the first row and column, add a row of —l's and 
then replace —1 by throughout. The result is BA(j),u — 1,2; ^jp). □ 

7. Comments. By exploiting properties of the dual, optimality theory 
for resolvable designs with two blocks per replicate has been developed. 
Section 3 establishes conditions for Schur-optimality of equal concurrence 
designs depending on the block discordancy coefficient 7. For 7 where Schur- 
optimality is not established, the class of competitors has been significantly 
narrowed via Schur-ordering for small r. Section 4 characterizes the class of 
all E-optimal designs whenever any group-affine design exists and further 
determines the Schur-best of the E-optimal designs. Sections 5 and 6 apply 
these results for the important cases k\ — k2 < 2, including explicit design 
constructions employing an equivalence with balanced arrays. It is evident 
from the constructions that many other designs, eliminated by the Schur- 
domination argument in Theorem 15, do exist. 

For large (r > v) replication, the problem can be quite different and op- 
timality work would proceed based on treatment, rather than block, con- 
currences. Here is a simple construction for that case: given a BIBD for v 
treatments in 6 blocks of size k, there is a Schur-optimal resolvable design in 
T>(v,b;k,v — k) formed by the blocks of the BIBD and their complements. 
While it is possible that these designs are equivalent to balanced arrays 
(when the starting BIBD is symmetric), it is also true that the resulting 
4>dhh' can be widely dispersed, for they are determined by the block concur- 
rence counts for the underlying BIBD. 
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The optimality route is clear for r < v — 1 and 7 < 7^ but f° r larger 7 
there are still unanswered questions. For instance, A- and E-optimality need 
not coincide and the problem of determining an A-best design remains open. 
For r < 4, the authors have solved the A-optimality problem in its entirety, 
including construction, and for larger r, have done this for the special cases 
of k\ — k 2 < 2. These and related results will be reported elsewhere. 



APPENDIX A: PROOFS 

1/2 

A.l. Proof of Lemma 1. Let D s be the diagonal matrix of square roots 

— 1/2 1 

of block sizes and write B d = N d D s . Multiplying C d by - and right and 

— 1/2 

left multiplying Cdual by D s , equations (1) and (2) become 

(16) \ C d = I~ - r N d D~ l N' d = 1- - r B d B' d = C d 
and 

(17) J D7 1 /2 CdualjD7 l/2 = j _ l _ D ~l/2 KNdD ~l/2 = j _ \ B , dBd = ^ 

Since (1) and (16) differ only by a constant and since the nonzero eigenvalues 
of B d B' d and B' d B d are identical, the eigenvalues of C d can be found from 
those of Cj ua | in (17), which depends on block concurrences and block sizes. 

Continuing, write C* dual = 1 - \A d for A d = B' d B d = DJ 1/2 N d N d D7 1/2 . 

1/2 1/2 

Regardless of the design d, A d D s 1 = rD s 1; r is an eigenvalue of A d 
corresponding to the zero eigenvalue common to Cdual and C d . One term of 
the spectral decomposition of A d is then 



(18) r(Dl /2 l)(Dl /2 iy _ 1 



(dI /2 i)'(dI /2 i) (h + k 2 ) 



J 



( ki Vhk 2 
k 2 



where J is an r x r matrix of ones. Subtract (18) from A d and denote the 
result A* d . Then a bit of manipulation, employing the fact that all four block 
concurrence counts for blocks in replicate h with blocks in replicate h! are 
determined by ^dhh' gives 

Since the eigenvalues of the 2x2 matrix in (19) are and k± + k 2 , the 
b = 2r eigenvalues of A* d are r copies of and times the r eigenvalues 
of M d . The eigenvalues of C d as stated in the lemma are now immediate. 
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A. 2. Proof of Lemma 17. The lower bound on e d \ follows from 
e<£L = max x'M d x 

x'x=l 



max x 

x'x=l 



> max x' 
x'x=l 
x'l=0 



P 



-+ 1 jI- 1 J + A d 
+ 7 ]I- 7 J+ A d 



P 



+ 7 + max x'A d x 

V x'x=l 
x'l=0 



P 



+ 7 + max x PA d Px 

V x'x=l 
x'l=0 



+ 7 + max x'PA d Px 

V x'x=l 

+ 7 + u dl . 



The penultimate equality holds since u d \ > and 1 is an eigenvector of 
PA d P with eigenvalue 0. Likewise, 



e dr = min x'M d x 

x'x=l 



mm x 

x'x=l 

p 



V 



-+7 7-7J + A d 



< — h 7 + min x P A d Px 

V x'x=l 



P 



-+J + u dr . 

V 

The last equality holds provided u dr < 0, for reasons similar to those above. 



If u dr > 0, the bound still holds, since e dr < 



tr(A/„) 



2 <r E 

V 



A. 3. Proof of Lemma 19. Suppose that A — J is of the suggested form. 

Write t(i) for ti + 1 2 H Hi and i (0 ) = 0. Then PAP = P(A- J)P, implying 

that 

maxeig(PAP) = max x' P{A — J)Px 

x'x=l 

= max y {A — J)y 

y=Px,x'x=l 



max 

y=Px,x 



t / %) \2 

!U-E( E vi) > 

i=i \j=t {i _ t) +i / 
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which is clearly nonpositive. 

Now, subscripting by the dimension, suppose that P r A r P r is nonpositive 
definite (npd). Exhaustive enumeration shows that A r — J r must have the 
form (13) for r = 3, 4, 5. Assuming nonpositivity implies that the form must 
hold for a given r, it will be shown that the same implication holds for r + 1. 
Denoting the upper left-hand side r x r submatrix of A r+ \ by A r (which 
is itself an adjacency matrix), it is claimed that P ri A T+ \P r+ i is npd =>• 
P r A r P r is npd. If not, then there exists x such that x' P r A r P r x > 0. Write 
y' = (x'P r ,0). Then y'l = => y' P r+ iA r+ iP r+ \y = y' A r+ \y = x'A r x > 0, a 
contradiction. 

Since P r A r P r is npd, the induction hypothesis says that 



. 1 i 



J, 



r+1 



-Jt, 




\ 


-Jt 2 










a 




— Jtn 




a' 


-1/ 



for some vector a rx \ of 0's and — l's. Indeed, by the induction hypothesis, 
every s x s principal minor of A r+ i — J r +i must have the form (13), so a 
may be partitioned as a' = (a^a^, ■ ■ ■ ,a' n ), where is either — l ti or 0^. If 



n 



1, the proof is done. If n = 2 and 02 = — , then 



A 



r + 1 



J, 



r+1 




a\ 

-Ju —h,< 



■V, 



>2 



1 




If now a\ = — lfj, then A r+ \ — J r +i contains a principal minor, 



Ma 



which is not of the form (13), contradicting the fact that the result holds 
for r = 3. Thus, at most one of 0,1,0,2 is nonzero. It follows immediately 
that for n > 2, the same statement holds for any pair Oi,Oi>. Permuting 
rows and columns as needed, it can be assumed that a, = 0^ for i < n and 
consequently that A r+ \ — J r +i has the form (13). 



A. 4. Eigenvalue equations for Hd and proof of Lemma 20. Any vector of 
the form (c'i, c' 2 , . ■ . , c' n ), where Cj € K'* is either a contrast vector or the zero 
vector, is an eigenvector of Hd with eigenvalue 0. Consequently, all other 
eigenvectors are of the form e' = {xil' t ,x<i\! t , . . . ,x n V t ) for some scalars 
X\, . . . , x n . The first equation in the system Hde = Xe is 7*1X1 — (1 — 7^2X2 — 
• • • — (1 — 7)t n x n = \x\; the other equations may be written similarly in 
order to see that Hde = Xe are equivalent to [D t — (1 — j)l n t']x = Ax, where 
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t = (ti, . . . , t n )' and x = (xi, . . . , x n )' . Thus, the remaining eigenvalues of Hd 
are the right eigenvalues of Df — (1 — 7)l n £'. Now, Dt is positive definite and 
E is nonnegative definite for 7 > so both have symmetric square root 
matrices and Dt is invertible. Thus, 

|A-(1-7)M'-A/|=0 

|A-(l-7)lni / -A/||A~ 1 |=0 

|#- ai^ -1 ! = o 

\D l t /2 ED] 12 - XI\ = (proving Lemma 20) 
|£ 1/2 A^ 1/2 - A/| = (needed in Theorem 15), 

the last step following because GG' and G'G have the same eigenvalues for 
any square G. Now, in the proof of Theorem 15 some elements of t are 
allowed to be zero (without loss of generality, t± = ■ ■ ■ = t z = for integer 
z > 1), in which case is not invertible. In this case there are z additional 
zero eigenvalues plus a reduced system of n — z equations in t z+ i, . . . ,t n . It 
is easy to see that \d\^ ' ED\^ 2 — AJ| = \Q z ,z — ^Iz\\Dl EiD\ ' 2 |, where E and 
Dt are the lower-right submatrices of E and Dt of order n — z. Thus, the n 
eigenvalues sought are still those of d\ I2 ED\ 12 and thus of E l / 2 D t E 1 l 2 . 



APPENDIX B: UNIFORM GBAS WITH FOUR GROUPS 



GBA(9,(1,1,1,2),2;2 



GBA(13,(1,1,1,2),2; 



3) 



11111 
00000 
01001 
01010 
00110 
10010 
00101 
11100 
10001 



11 1 11 
11111 

01100 
0001 
10010 
0000 1 
1000 1 
00110 
10100 
00101 
01010 
1001 
11000 
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GBA(17, (2, 2, 2, 3), 2; 4) GBA(21, (2, 2, 2, 3), 2; 5) 



111111111 


111111111 


000000000 


111111111 


100100101 


110000110 


001101100 


101010100 


100011001 


011010001 


101110010 


000101100 


010010110 


001010010 


010111100 


100000011 


100001110 


010101010 


111001001 


010000101 


101010100 


100101001 


000101011 


011001001 


011001010 


101001001 


110100010 


010011010 


010110001 


011100100 


001010011 


001001110 


011000101 


101100010 




100011100 




000110101 




110110000 




000110011 
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